Scaling of the superfluid density in superfluid films 
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We study scaling of the superfluid density with respect to the film thickness by simulating the x — y 
model on films of size L x L x H {L » H) using the cluster Monte Carlo. While periodic boundary 
conditions where used in the planar (L) directions, Dirichlet boundary conditions where used along 
the film thickness. We find that our results can be scaled on a universal curve by introducing an 
effective thickness. In the limit of large H our scaling relations reduce to the conventional scaling 
forms. Using the same idea we find scaling in the experimental results using the same value of 
V = 0.6705. 
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Scaling is a central idea in critical phenomena near a 
second order phase transition and in field theory when we 
are interested in the continuum limit In both cases 
we are looking at the singular behavior emerging from 
the overwhelmingly large number of degrees of freedom, 
corresponding to the original cutoff scale, which need to 
be integrated out leaving behind long-wavelength degrees 
of freedom which vary smoothly. Their behavior is con- 
trolled by a dynamically generated length scale, the cor- 
relation length ^. Such a fundamental idea is difficult to 
test theoretically because it requires a study of an over- 
whelmingly large number of interacting degrees of free- 
dom. Experimentally, however, one hopes to be able to 
study scaling in finite-size real systems near a second or- 
der phase transition. Namely, the system is confined in 
a finite geometry (for example, film geometry) and the 
finite-size scaling theory is expected to describe the be- 
havior of the system near the bulk critical temperature 
T\. Liquid '^He has been a good real system for test- 
ing finite-size scaling theory and measuring the critical 
exponents that are associated with the most singular be- 
havior in thermodynamic quantities near T\. However, 
measurements of the superfiuid density and the spe- 
cific heat 1^ on helium films fail to verify the finite-size 
scaling theory. 

The situation of the specific heat has been recently 
clarified where the choice of the boundary condi- 

tions was a key factor in comparing the more recent mea- 
surements of the universal function |Q and that obtained 
theoretically. On the other hand, while new experiments 
for the specific heat under confined geometries have been 
planned to be conducted under more ideal microgravity 
conditions § , the problems related to the measurements 
of the superfluid density |^ are still outstanding. 

In this paper we use the x — y model and the cluster 
Monte Carlo method to calculate the superfluid density 
on films of size LxLxH {L » H) with periodic bound- 
ary conditions in the planar L-directions and Dirichlet 
boundary conditions (vanishing order parameter) along 
the film thickness dimension. The same model, geome- 



try, and boundary conditions where used in Ref. |^] to 
calculate the specific heat, and a very good agreement 
between the theoretically calculated and experimentally 
determined universal functions was found. In this paper 
we show that the superfluid density is a far more sen- 
sitive observable than the specific heat with respect to 
the requirement that one needs to use very thick films 
{H oo) to verify scaling with respect to the film thick- 
ness H . We have found that in order to achieve scaling 
for rather small values of H (as in the case of the spe- 
cific heat) we need to modify the scaling expressions by 
using a concept of an effective thickness. Our introduc- 
tion of an effective thickness i?e// = H + D (where D 
is a finite dynamically generated length scale) is no vi- 
olation of scaling, since at large Heff the constant D 
can be neglected. Scaling for all the values of H used 
in our calculation is achieved with the expected value of 
ly = 0.6705 |9). A similar modification to the scaling for- 
mula allows scaling of the experimental results of Rhee 
et al. 1^ for the superfluid density with the same value 
of ly. 

For the x — y model on a lattice, the helicity modu- 
lus T^(T)/J as defined in Refs. is calculated as the 
ensemble average of j) cos{9i — Oj){e^ ■ e^)^ — 

f}Cl2{ij) sm{6i — Oj)efi ■ Cij)^). Here V is the volume of 
the lattice, (3 = J/ksT, is the unit vector in the corre- 
sponding bond direction, and eij is the vector connecting 
the lattice sites i and j. In the following we omit the vec- 
tor index since we will always refer to the x-component of 
the helicity modulus and due to the isotropy T^, = Ty. 
The connection between the helicity modulus and the 
superfluid density ps is established by the relation |l^] 
Ps{T) = (m/?i)^T(r) where m denotes the mass of the 
helium atom. 

In Ref. we studied the helicity modulus T for the 
x — y model in a film geometry with periodic bound- 
ary conditions in the _ff -direction. In a certain tem- 
perature range around the bulk critical temperature Tx 
where the bulk correlation length £,{T) becomes of the or- 
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der of the film thickness H the quantity TH/T exhibits 
effectively two-dimensional behavior and a Kosterlitz- 
Thouless phase transition takes place at a temperature 
T^^{H) < Tx. We found that the critical temperature 
T^^iH) approaches Tx in the limit H 



T'^°{H) = Tx[l 



(1) 



where the critical exponent v is the same as the experi- 
mental value V = 0.6705 g and the value Tx/ J = 2.2017 
I p^ . We also demonstrated that TH/T is a function of 
the ratio H/£^{T), i.e. the dimensionless quantity 



T{T,H)H 
T 



<i>{tH 



(2) 



is a function of x = tH^/'^ only. We found that when 
we plotted the calculated T{T, H)H/T as a function of 
X, in the limit L — + oo, our results for all thicknesses 
H collapse on the same universal curve. Thus, simple 
scaling holds for periodic boundary conditions. 

In this paper we consider periodic boudary conditions 
in the planar L-directions and Dirichlet boundary con- 
ditions along the thickness direction. Fig.|l| displays our 
Monte Carlo data for the helicity modulus in units of the 
lattice spacing a and the energy scale J for the film of 
fixed thickness H — 4. Dirichlet boundary conditions 
strongly suppress the values of the helicity modulus as 
compared to the case of periodic boundary conditions 
along H. As a consequence, films with Dirichlet bound- 
ary conditions have lower critical temperatures than films 
with periodic boundary conditions. 




FIG. 1. The helicity modulus T(r, H, L) as a function of T 
for various lattices x 4 with Dirichlet boundary conditions 
(Dbc) in the ii"-direction. 

We eliminate the finite size effects in the L-direction by 
studying the quantity K = T / (TH) in the L oo limit. 



For a fixed thickness H and at temperatures T below 
but sufficiently close to the critical temperature T^^ {H) , 
the system behaves effectively two-dimensionally [|l2|,0 . 
It was demonstrated in Ref. that we can use the 
Kosterlitz-Thouless-Nelson renormalization group equa- 
tions | p5| to derive an expression for the planar L- 
dependence of K. This expression can be used to ex- 
trapolate the computed values K{T, H, L) obtained on 
lattices of finite L to the L = oo limit, for a fixed H. 
This has been clearly demonstrated in Ref. for the 
case of periodic boundary conditions. The quality of our 
extrapolation is the same as in Ref. and we omit such 
demonstration here due to lack of space. In the following 
we shall drop the dependence of T on L implying that 
we refer to the extrapolated L ^ oo values. 

In Fig.| we plot T{T,H)H/T versus tH^/'' for the 
thicknesses H = 12, 16, 20, 24 to check the validity of 
the scaling form (H) using the experimental value of 
h' — 0.6705 |9|. We do not obtain a universal scahng 
curve, thus scaling according to the expression is not 
valid for the films with thicknesses up to iJ = 24. 



2.0 



X 



1.0 



0.0 



7 H=12 
H=16 
□ H=20 
c H=24 



-6.0 



-5.0 



-4.0 
tH" 



-3.0 



-2.0 



FIG. 2. r{T,H)H/T as a function of ti/^/" for various 
thicknesses, u — 0.6705. 

Let us, therefore, pursue another line of thought. 
In Fig.^ we show the layered helicity modulus 
Ti(2;)/J,where z counts the layers, computed on a 60 x 
60 X 20 lattice at the temperature T/J = 2.1331. The 
quantity T l{z) / J is just the helicity modulus determined 
for each layer separately. The layered helicity modulus 
is symmetric with respect to the middle layer where it 
reaches its maximum and decreases when the bound- 
aries are approached. Although the helicity modulus 
T(T, H,L)/ J is not the average of the quantity (z) / J 
over all layers, the curve in Fig. ^ is an approximation 
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to the profile that the superfluid density develops in thin 
films. The basis for the standard scaling argument is the 
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FIG. 3. The approximate profile Tl{z) of the helicity 
modulus computed on a 60 x 60 x 20 lattice at T = 2.1331, 
i.e. close to the critical temperature, T^^{20) — 2.1346. 



following. For large H and very close to the critical point 
where £,{T) is very large, the "penetration" depth A(T) 
of the superfluid density inside the film is of the order of 
the correlation length. Thus, in the limit where all other 
length scales are small compared to H and ^, if we plot 
YH versus z/H (or z/£,) we should find scaling. However, 
for small H there is at least one length scale D (which 
for H ^ D needs to taken into account) which has the 
following origin. The length scale D contains information 
on how fast T(z) rises from T(z = 0) = 0. Namely the z- 
derivative of T{z) is not universal, it depends on how we 
have imposed the Dirichlet boundary conditions. There 
are many ways to make the order parameter vanish at 
the boundary. It can be made to be zero when averaged 
over a boundary area A — I x I. In our case of staggered 
boundary conditions I = \/2. If we had chosen Dirichlet 
boundary conditions where the order parameter is zero 
over an area with I > \f2 we would have found a slower 
rise of T(z) from its zero value at the boundary. If this 
initial "faster rise" of the superfluid density is neglected, 
the rest of T(z) can be fit to Acosh(z/^) + B with only 
one length scale, the correlation length ^. Thus the curve 
y(z) can be thought of as made of two contributions, and 
scaling at small values of H can be obtained only if the 
film size is extended. We can imagine that this thinner 
film of size H is obtained from a thicker one by a process 
of forcing the superfiuid density to go to zero faster than 
its "natural way" by a "speed" dictated by the severity 
of the boundary conditions. 



The lack of scaling with the expected critical expo- 
nent V = 0.6705 indicates that the critical temperatures 
Tl^{H) do not satisfy Eq. (0). Because of the argument 
given earlier about the profile of the superfluid density 
we may expect an effective film thickness i?e// to en- 
ter the scaling expressions (0) and (H). The simplest 
assumption is i?e// = H + D where D is a constant. In- 
deed by replacing H with H^f / in these equations for the 
film thicknesses H = 12, 16, 20 we obtain Xc = —3.81(14) 
and D = 5.79(50) with v = 0.6705. In Fig.| we plot 
T{T, H)Heff /T as a function of tH^j:'^ for films with 
H = 12, 16, 20, 24 where v = 0.6705. The data for the 
helicity modulus collapse onto one universal curve. 
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FIG. 4. T(r, H)H^ff/T as a function of tHl^/j^ 
thicknesses. H^ff — H + 5.79 and u — 0.6705. 
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We wish to test the assumption that the boundaries in- 
troduce an effective thickness into the scaling expression 
(^ further. Janke and Nather studied the thickness 
dependence of the Kosterlitz-Thouless transition tem- 
perature of the Villain model with open boundary condi- 
tions (interactions of the top and bottom layer only with 
the interior film layers). They found, however, that in or- 
der for scaling to occur they needed to use a value for ly 
higher than the value believed for the model. We replace 
H in Eq. (|l|) by the effective thickness i?e// — H + Dy. 
Indeed, taking the expected value u = 0.6705 we find 
Dv = 1.05(2) and Xc = —1.62(2) and a good quality of 
fit. We can understand the increment D as an effective 
scaling correction which renders the scaling relations (|l|) 
and (^ valid even for very thin films. For large thick- 
nesses H the increment D can be neglected and we re- 
cover the conventional scaling forms. This result means 
that the film thicknesses considered in Ref. [|l^ were still 
too small to extract the expected value of the critical 
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exponent v from the iJ-dependence of the critical tem- 
perature (|^) without the help of an effective thickness 

In the experimental situation it is possible to imagine a 
similar situation where a length scale D emerges and cor- 
responds to an average defect distance on the substrate. 
In Fig.H we achieve approximate collapse of the data for 
the superfluid density ps for films of various thickness d 
(d is in \irn) given in Refs. by plotting ps{t,d)def f / p 
versus ^cJg^j with v = 0.6705 and deff = d + 0.145. We 
obtained the effective thickness by examining the reduced 
temperatures tfg{H) where finite-size effects set in. Ac- 
cording to finite-size scaling theory tfg has to fulfill the 

relation tfs oc d~^^'^, thus in our case tfs oc t^e//"- T^^^ 
data points corresponding to the film with d — 3.9fim 
deviate from the universal curve; we attribute this to the 
anomalous behavior of these data. Namely, in general 
\tfs{di)\ > \tfs{d2)\ if di < d2, but this is not the case 
for c?i = 2.8/im and o?2 = 3.9/im (cf. Refs. [^). The scal- 




FIG. 5. Scaling of the superfluid density data of Rhee et 
al. with the effective thickness deff — d+0.145. u = 0.6705 
and all lengths are in pm. 

ing of the experimental data with the aid of this finite 
length scale indicate that the average distance between 
defects will introduce a new length scale from which one 
has to stay away in order to see simple finite size scal- 
ing. Theoretically there are two different limits in which 
simple scaling with H should be expected to be valid. 
The clean boundary limit, ^ ^ H D, and the dirty 
boundary limit where Z? <C ^ ~ i?. 

In conclusion, we find that the results of our simula- 
tions for the superfluid density on rather small thickness 
films obey a scaling relation where an effective thickness 
is introduced. Conventional scaling relations is expected 



to be valid on films of much larger thickness. Applying 
the same idea of the effective thickness on the experi- 
mental data, we found that the long standing problem of 
lack of conventional scaling in the data of Rhee et al ||] 
can be resolved in a simple way without resorting to any 
departure from scaling nor to using unrealistic values for 
v. Clearly, more experiments with different substrates 
are desirable. 
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